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Abstract 


A negative correlation was found to exist between the low-energy spectral index and the redshift of gamma-ray 
bursts (GRBs) by Amati et al. It was later confirmed by Geng & Huang and Gruber et al., but the correlation was 
also found to be quite dispersive when the sample size was significantly expanded. In this study, we have 
established two even larger samples of GRBs to further examine the correlation. One of our samples consists of 
316 GRBs detected by the Swift satellite, and the other one consists of 80 GRBs detected by the Fermi satellite. It 
is found that there is no correlation between the two parameters for the Swift sample, but there does exist a weak 
negative correlation for the Fermi sample. The correlation becomes even more significant when the spectral index 
at the peak flux is considered. It is argued that the absence of the correlation in the Swift sample may be due to the 
fact that Swift has a very narrow energy response so that it could not measure the low-energy spectral index 
accurately enough. Further studies based on even larger GRB samples are solicited. 
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1. Introduction 


Gamma-ray bursts (GRBs) are the most violent stellar 
explosions. They were initially discovered in 1967 (Klebesadel 
et al. 1973) by the Vela Satellites. Up to now, they have been 
intensively studied for more than forty years since more and 
more GRBs are discovered (Qin et al. 2021; Yuan et al. 2022a; 
Liu et al. 2022). Generally, GRBs are believed to be triggered 
by the death of massive stars or by the merger of binary 
compact stars. The empirical correlation between different 
parameters of GRBs is one of the present interesting objects. 
For instance, the correlation between E, (the peak energy) and 
Eiso (the equivalent isotropic energy), i.e., the so-called Amati 
relation is investigated by Amati et al. (2002, 2009), Virgili 
et al. (2012), Azzam & Alothman (2013), Geng & Huang 
(2013) and Demianski et al. (2017). The correlation between 
E, (the peak energy) and L, (the peak luminosity), i.e., the so- 
called Yonetoku relation is investigated by Yonetoku et al. 
(2004), Ghirlanda et al. (2005) and Zhang et al. (2012). 
Yonetoku et al. (2010) and Tsutsui et al. (2013) have discussed 
both relations. Furthermore, Wang et al. (2020) has even 
discussed other more empirical correlations beyond the above 
relations. 

The redshift (z) is an important parameter of GRBs. Many 
people have studied the relation between z and other 
parameters. For example, Zhang et al. (2014a) investigated 
the correlations between z and Liso, Eiso Ep,est. Ep,obs 


(Ep,obs = Ep,rest/(1 + Z)) (also see Sakamoto et al. 2011; 
Ukwatta et al. 2012). Especially, they found a positive 
correlation between z and Liso for the Swift GRBs. Wei & 
Gao (2003) and Zitouni et al. (2014) studied the z-E;,, relation 
and the z-L;,. relation. Moreover, the relation between z and L, 
has also been investigated by many researchers (Lloyd- 
Ronning et al. 2002; Goldstein 2012; Salvaterra et al. 2012; 
Zhang et al. 2014b; Zitouni et al. 2018). 

The low-energy spectral index, a, is also an important 
parameter that characterizes the prompt emission of GRBs. 
Geng et al. (2018) have tried to derive the value of a through 
numerical simulations by considering the synchrotron emission 
mechanism. Especially, the correlation between a and E, has 
been studied by many researchers (Gruber et al. 2014; Li et al. 
2019; Duan & Wang 2020; Li 2022). Tang et al. (2019) also 
investigated the correlation between a and Eiso- 

Interestingly, Amati et al. (2002) found that there seems to 
be a negative correlation between z and a, which reads 
loga = (—0.78 + 0.13)log(1 + z) + (0.39 + 0.04). Here the 
symbol log means the logarithm base is 10. Their study is 
conducted on a sample that only includes nine GRBs. Later, 
Geng & Huang (2013) and Gruber et al. (2014) confirmed the 
existence of such a correlation, which was updated as 
loga = (—0.42 + 0.07) log(1 + z) + (0.11 + 0.02). Their 
sample consists 65 GRBs, which is much larger than that of 
Amati et al. (2002). However, Geng & Huang (2013) also 
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noted that the correlation is quite dispersive. It is thus necessary 
to further examine the existence of such a z — a correlation. In 
this study, we will further expand the GRB sample to include 
316 GRBs detected by Swift and 80 GRBs detected by Fermi. 
We then use this significantly expanded sample to check the 
correlation between z and a. 

The structure of this article is organized as follows. A 
detailed description of the selection of the data sample is 
provided in Section 2. The correlation between the low-energy 
spectral index and the redshift, and some other relations 
between various parameter pairs, are explored in Section 3. 
Finally, our conclusions and discussion are presented in 
Section 4. 


2. Sample 


GRBs detected by Swift and Fermi are used in this study. 
Two conditions are applied in selecting the appropriate GRBs. 
First, the redshift of the burst should be available. Second, the 
spectrum should be well defined. The time-averaged spectra of 
GRBs are usually fitted with three kinds of functions: a single 
power-law function, a cutoff power-law function, and the so- 
called Band Function (Band et al. 1993). Most of the GRBs 
spectra can be well fitted by the Band function, which is 
expressed as 


F(E) 
E \@ Q- a)E —(a— B)Ep | 
A Ga) exp |- E, |. E< Q-a ° 
E E \-8 —@-8)E, CEO —(a — B)Ep 
(mow) exp(a — o| 100 keV(2 — | a 
(1) 


where E is the photon energy. There are four parameters in this 
equation: the scaling factor (A), the low- and high-energy 
spectral indices (a and (3, respectively), and the peak photon 
energy (£,). In our notation, œ and 8 are positive by using their 
absolute values. We have collected all the spectrum parameters 
of those Swift and Fermi GRBs with the redshift measured. 

Fermi has a very wide energy band, i.e., 8—35,000 keV. So, 
Fermi GRBs are usually best fitted by the Band function, which 
is represented by the lower-energy spectral index (aœ) and the 
high-energy spectral index (3). On the other hand, the Swift/ 
BAT detector has a very narrow energy response of 
15-150 keV so that Swift GRBs are generally best fitted by a 
single power-law function or by a cutoff power-law function. 
In these cases, the derived power-law index could be regarded 
as a useful representation of the low-energy spectral index (a), 
since it is measured in the soft y-ray range. Note that for the 
Swift GRBs, the 3 parameter is completely unavailable. 

For the Swift GRBs, the relevant data are acquired by 
inquiring on the NASA Swift website. As a result, a total 
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Figure 1. Swift sample on the log a—log(1 + z) plane. Here a is the low- 
energy spectral index and z is the redshift. The solid line is the best fit result 
and the shaded region shows the corresponding lo range. 


number of 316 GRBs are included in our Swift sample, all of 
which have the necessary redshift and spectrum data. The time 
span of these GRBs ranges from 2005 January 26, to 2023 
January 16. 

For the Fermi GRBs, the data are collected mainly by 
consulting NASA’s HEASARC (the High Energy Astrophysics 
Science Archive Research Center) database. As a useful 
supplement, 24 bursts were taken from Jochen Greiner’s online 
GRB catalog.’ Finally, our Fermi sample consists of 80 GRBs. 
The time span of these Fermi GRBs ranges from 2008 
September 5 to 2018 July 20. 


3. Correlations 


In this section, the Swift and Fermi GRB samples are 
employed to explore the empirical correlations concerning a 
and z. Figure 1 shows the Swift sample on the 
log a — log(1 + z) plane. Generally, we see that the data points 
are quite scattered and there is no obvious correlation between the 
two parameters. The solid line in Figure 1 is the best-fit result of 
the data points, which corresponds to loga = (0.029 + 0.030) log 
d+z)+(0.15+0.014). But note that it has a very small 
correlation coefficient of r=0.055, which indicates that essen- 
tially no correlation exists. Here, the error ranges of the best-fit 
parameters are given in lo ranges throughout this paper. 

In Figure 2, we plot the Fermi GRB sample on the 
log a—log(1 + z) plane in the left panel. The observational data 
points are also fitted with a linear function. Note that GRB 
120712A seems to be an obvious outlier on the plane, whose 
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Figure 2. The spectral index plotted against redshift for the Fermi GRB sample. The left panel shows the low-energy spectral index (œ) vs. redshift (z). The right panel 
shows the high-energy spectral index (8) vs. redshift (z). The solid lines correspond to the best fit results to the data points and the shaded regions show the 
corresponding lo ranges. Note that in the left panel, the isolated data point (GRB 120712A) on the top is not included in the linear fit since it is an obvious outlier. 
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Figure 3. a vs. E, (left panel) and E, vs. z (right panel) for the Fermi GRBs. The solid lines correspond to the best fit results to the data points and the shaded regions 
show the corresponding lo ranges. Note that in the left panel, the isolated data point (GRB 120712A) on the top is not included in the linear fit since it is an obvious 


outlier. 


low-energy spectral index is a=6.1+0.67. We exclude this 
event in the fitting procedure. Interestingly, we could see that there 
is a negative correlation between a and z. The best-fit result is 
loga = (—0.18 + 0.11)log(1 + z) + (0.012 + 0.049), with 
a correlation coefficient of r= -— 0.19. But note that this 
expression, especially the slope of the line (i.e., —0.18), is quite 
different from those of previous studies. Also, the correlation here 
is generally very weak. As a comparison, the right panel of 
Figure 2 shows the Fermi GRB sample on the log G—log(1 + z) 
plane. We see that there is no correlation between the two 
parameters, which is consistent with previous results of 
Geng & Huang (2013). The best-fit line corresponds to log3 = 
(0.015 + 0.055) log(1 + z) + (0.36 + 0.025), with a very small 
correlation coefficient of r= 0.031. Comparing the two panels of 
Figure 2, it is quite clear that the low-energy spectral index (œ) is 


still obviously different from the high energy index (3) when they 
are plotted against the redshift. 

The Fermi GRBs have well-measured E, data. So, in 
Figure 3, we plot them on the a — E, plane and the E, — z 
plane. We see that there is no obvious correlation either 
between a and E,, or between E, and z. The results are also 
somewhat different from those of previous studies. For 
example, Geng & Huang (2013) argued that E, and z are 
positively correlated. Note that GRB 120712A is again an 
obvious outlier in this figure and is not included in the fitting 
procedure. 

Since the peak time of the flux is an important stage of a 
GRB, we have also investigated the features of the peak flux 
parameters. For this purpose, we could only adopt the Fermi 
GRB sample, since peak flux parameters are not available for 
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Figure 4. The peak flux spectral indices plotted against the redshift for the Fermi GRB sample. The left panel shows the low-energy spectral index at the peak flux 
moment (Qpeak) VS. the redshift. The right panel shows the corresponding high-energy spectral index (Gpeax) vs. the redshift. The solid lines correspond to the best fit 
results to the data points and the shaded regions show the corresponding lo ranges. 


many Swift GRBs. Here, the low-energy spectral index at the 
peak flux of the burst is denoted as Opeax. Correspondingly, the 
high-energy spectral index and peak photon energy at the peak 
flux are denoted as peak and Epeak» respectively. 

In the left panel of Figure 4, the Fermi GRB sample are 
plotted on the log dear — log(1 + z) plane. Here, we could 
see that there is an obvious negative correlation between Qpeax 
and z, although it is still quite scattered. The best-fit result is 
log Qpeak = (—0.29 + 0.15)log(1 + z) + (—0.042 + 0.059), 
with a correlation coefficient of r= — 0.29. Comparing with 
the slope of —0.42 previously derived for the loga — 
log(1 + z) correlation (Geng & Huang 2013), here the 
slope of the log Qpeak—log(1 + z) relation (i.e., —0.29) is 
significantly flatter. The right panel of Figure 4 shows the 
Fermi GRB sample on the log peak — log(1 + z) plane. There 
is also a weak correlation between the two parameters, 
which can be best fitted as log Geax = (—0.088 + 0.070) log 
(1 +z)+ (0.39 + 0.027), with a correlation coefficient of 
r=-—0.20. Note that the data points in the right panel are 
much dispersive as compared with those in the left panel. 

Figures 2 and 4 show that there is a weak correlation 
between the low-energy spectral index and the redshift. 
Especially, the a—z relation is much less significant than the 
Qpeak-Z relation. A natural speculation is that a and Qpeax 
should be positively connected, then the above two relations 
should be largely similar. The reason that leads to such a 
difference thus deserves to be examined. Figure 5 plots a 
against Qpeak for the Fermi sample. The solid line shows the 
case when a equals Qpeax. We see that a and Qpeak are not 
strictly connected, which could explain their different depend- 
ence on the redshift. It reflects the fact that the y-ray spectrum 
is highly variable during a GRB. 

Figure 6 illustrates the relations between Qpeax, Epeak and 
z. There is a weak correlation between Epeak and z, which reads 


0.4 0.8 1.2 1.6 
Apeak 


Figure 5. aca, plotted vs. a for the Fermi GRB sample. The solid line 
corresponds to the case of a = peak. 


log Epeak = (0.26 + 0.33)log(1 + z) + (2.30 + 0.13), with 
the correlation coefficient being r= 0.13. Similarly, there is also 
a weak correlation between peak and Epeak, Which corresponds 
to log Qpeak = (—0.072 + 0.071) log Epeak + (0.0088 + 0.17), 
with the correlation coefficient being r= — 0.16. 

Swift/BAT has a relatively narrow passband, i.e., 
15-150 keV. Therefore, the measured spectral index should 
correspond to the low-energy spectral index of the Band 
function. Some GRBs are simultaneously detected by Swift and 
Fermi satellites. It is then interesting to know whether the 
spectral index measured by Swift is consistent with that 
measured by Fermi. In Figure 7, we have screened out all the 
overlapping GRBs between the Swift sample and Fermi 
sample, and compared their spectral indices. Here, œ” is the 
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Figure 7. The low-energy spectral index measured by Fermi (a^), plotted vs. 
the spectral index measured by Swift (a S), Note that only overlapping GRBs 


between the Fermi sample and the Swift sample are shown here. The solid line 


corresponds to the case of a” = añ. 


spectral index measured by Swift, while a” is the low-energy 
spectral index measured by Fermi. We see that aŭ and a” are 
not equal for each event. The former is systematically larger 
than the latter. Also, the data points are quite scattered. For 
these overlapping GRBs, we have also compared their a% with 
their low-energy spectral index at the peak flux as measured by 
Fermi (Qheak): The results are shown in Figure 8. Similarly, we 


see that a“ again is generally larger than abeak Figures 7 and 8 
clearly show that different detectors could generate very 
different results for the spectrum of even the same GRB, 
which indicates that acquiring the spectra of GRBs is still a 
very difficult task. The different a — z correlations of different 
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Epeak (left panel), and Epeak plotted vs. z (right panel) for the Fermi GRB sample. The solid lines correspond to the best fit results and the 
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Figure 8. The low-energy spectral index at the peak flux measured by Fermi 
(kea) plotted vs. the spectral index measured by Swift (añ). Note that only 


overlapping GRBs between the Fermi sample and the Swift sample are shown 


here. The solid line corresponds to the case of OFF aie =a’, 


GRB samples could thus be caused by the systematic distortion 
in measuring the spectral indices. 


4. Discussion and Conclusions 


In this study, we use the Swift GRBs and Fermi GRBs to 
explore the possible correlation between the low-energy 
spectral index and the redshift. For the Swift GRB sample, it 
is found that there is no correlation between a and z (see 
Figure 1). On the other hand, there is a weak correlation 
between a and z for the Fermi GRB sample (see Figure 2). The 
correlation is even more obvious when the peak flux spectrum 
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is considered, i.e., when Qpeak is plotted versus z (Figure 4). 
The different features of the two samples may be caused by the 
different energy responses of the two detectors. The energy 
band of Fermi is very wide (8—35,000 keV), which ensures that 
it can present a much better description of the spectra of the 
detected GRBs. However, the Swift/BAT has a very narrow 
energy response (15-150 keV). As a result, the spectra of Swift 
GRBs are usually best fitted by a single power-law function or 
by a cutoff power-law function. In these cases, the power-law 
spectral index may significantly deviate from the true low- 
energy spectral index. This conjecture was confirmed when the 
spectra of the overlapping GRBs of the two samples were 
scrutinized. It is found that when a GRB is simultaneously 
detected by both Fermi and Swift, then the spectral index 
reported by Swift is usually quite different from the low-energy 
spectral index (a) reported by Fermi (Figures 7 and 8). It 
reflects the difficulty in spectral observations of GRBs. 

Amati et al. (2002) argued that the correlation between a and z 
is loga = (—0.78 + 0.13)log(1 + z) + (0.39 + 0.04). Note 
that their sample size is very small. Later, Geng & Huang (2013) 
updated the correlation as log a = (—0.42 + 0.07)log(1 + z)+ 
(0.11 + 0.02) using a significantly expanded sample. We see 
that the slope between loga and log(1 + z) becomes much 
flatter, and the correlation is also very dispersive. According to our 
current study, the correlation is now loga = (—0.18 + 0.11) 
log(1 + z) + (0.012 + 0.049) for the Fermi sample, and it is 
log Apeak = (—0.29 + 0.15)log(1 + z) + (—0.042 + 0.059) for 
the peak flux spectrum. The correlation becomes even flatter and 
more dispersive. Anyway, based on this study, we can conclude 
that a weak correlation does exist between a and z. Also, further 
study on this issue still deserves to be conducted in the future. For 
this purpose, detectors with a relatively wide energy response are 
necessary to accurately measure the exact spectra of GRBs. We 
note that the Einstein Probe, a wide field (3600 square degrees) 
X-ray satellite that is scheduled to be launched in late 2023 (Yuan 
et al. 2018, 2022b), might be very helpful in this aspect. Einstein 
Probe is very sensitive in X-rays and may markedly help increase 
the sample size of well-localized GRBs. 
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